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Abstract 

We study the large order behaviour in perturbation theory of the Bjorken, Elhs-Jaffe 
and Gross-Llewellyn-Smith sum rules. In particular, we consider their first infrared 
renormalons, for which we obtain their analytic structure with logarithmic accuracy 
and also an approximate determination of their normalization constant. Estimates of 
higher order terms of the perturbative series are given. The Renormalon subtracted 
scheme is worked out for these observables and compared with experimental data. 
Overall, good agreement with experiment is found. This allows us to obtain do and 
some higher-twist non-perturbative constants from experiment: do = 0.141 it 0.089; 
/3,Rs(l GeV) = -0.124l0il^ GeV^. 
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1 Introduction 



Deep Inelastic Scattering (DIS) is one of the few places where one can test, on a very 
solid theoretical ground, asymptotic freedom and the operator product expansion against 
experiment. This is so because, for the moments of the different structure functions, the 
transfer momentum lies in the euclidean region and far away of the physical cuts (for 
large). Therefore, the theoretical predictions do not rely on any kind of quark-hadron duality 
and may provide with solid and very clean determinations of and some non-perturbative 
matrix elements. Equally interesting is the study of the interplay between the perturbative 
and non-perturbative regime. This can only be done if a full control on the perturbative 
series is achieved. 

Within DIS, special consideration deserve the sum rules for which the matrix elements 
are related to symmetry generators (in this case they can be computed in absolute value 
within perturbation theory), like the Gross-Llewellyn-Smith (GLS) sum rule [1], or to some 
low energy constants (which can be directly measured from experiment), like the Bjorken 
sum rule [2]. In this paper, we will concentrate on the Bjorken, Ellis- Jaffe [3] and GLS 
sum rules. Their leading-twist term has been computed with next-to-next-to-leading order 
(NNLO) accuracy and they have been measured with increasingly good accuracy over the 
years. 

By the Operator Product Expansion, the short- and long-distance contributions are sep- 
arated and a sum rule, M, can be expressed in the following way: 

M = C<J> + S^ + ..., (1) 

with short-distance Wilson coefficients C, -B, ... and long-distance matrix elements <J>, 

<R>, The perturbative series of C are expected to be asymptotic and, therefore, to 

diverge for a high enough order in perturbation theory. Moreover, in schemes without strict 
separation of large and small momenta, such as MS, they are believed to be non-Borel 
summable. This is because when calculating the matching coefficients C, . . . , the integrals 
run over all loop momenta, including small ones. Therefore, they also contain, in addition to 
the main short-distance contributions, contributions from large distances, where perturba- 
tion theory is ill-defined. These contributions produce infrared renormalon singularities [4], 
factorially growing contributions to coefficients of the perturbative series, which lead to am- 
biguities ~ (Aqcd/Q)^" in the matching coefficients C, .... Similarly, matrix elements of 
higher- dimensional operators <R>, . . . also contain, in addition to the main large-distance 
contributions, contributions from short distances, which produce ultraviolet-renormalon sin- 
gularities. They lead to ambiguities of the order Aq^jj times lower-dimensional matrix 
elements (e.g., < J>). These two kinds of renormalon ambiguities should cancel in physical 
observables [5-10], in this case M. 

The intrinsic (minimal) error associated to the perturbative series is of the order of the 
higher twist correction. Thus, one can not unambiguously determine the higher twist terms, 
unless a prescription to deal with the perturbative series that has power-like accuracy is 
given. In this paper we will adapt to this case the prescription used for heavy quark physics 
in Refs. [11-13]. The idea is that the leading divergent behaviour of the perturbative series is 
related to the closest singularities in the Borel plane of its Borel transform. In heavy quark 
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physics, they he on the positive axis (infrared renormalons) . In the case of the sum rules 
considered in this paper, the closest singularities lie on the positive and negative axis at equal 
distance to the origin. We will assume that the one in the positive axis will dominate the 
asymptotics of the perturbative series. Since these singularities cancel against the ultraviolet 
renormalons of the low energy dynamics of the twist-4 operators, the proposal will be to shift 
the singularities from the perturbative series to the twist-4 operators. We will refer to this 
prescription as the Renormalon Subtracted (RS) scheme and apply it to the Bjorken, EUis- 
Jaffe and GLS sum rules. We will obtain the contribution of the leading infrared renormalon, 
subtract it from the perturbative series and add it to the low energy matrix elements, in our 
case, to the twist-4 operators. In this way one enlarges the range of convergence of the Borel 
transform of the perturbative series, which can be defined with power-like accuracy. This 
procedure has proven to be extremely successful in heavy quark physics, where it has been 
shown that perturbation theory works very well and good determinations of non-perturbative 
subleading corrections have been obtained using either lattice or experimental data. One 
should be aware, however, that in heavy quark physics one was in an optimal situation, since 
the singularity in the Borel plane was quite close to the origin {u = 1/2). We are now going 
to be in a less optimal situation, since the closest singularities lie at w = ±1. The physical 
situation is also completely different since now we are talking of a system made with light 
fermions. Therefore, it is interesting to investigate if a similar improvement is obtained in 
this case. We will do so in this paper. 

The paper is organized as follows. In the next section we will introduce the relevant sum 
rules. In Sect. 3 the Borel transform of the first infrared renormalon of the leading-twist 
Wilson coefficient will be calculated with leading log accuracy, as well as the normalization 
constant and estimates of the higher order terms of the perturbative series. The RS scheme 
will be worked out in Sect. 4. In Sect. 5 the comparison with the experimental data will 
be done allowing us the extraction of some non-perturbative matrix elements. Finally, the 
conclusions are presented in Sect. 6. 

2 Sum Rules 

The Bjorken and Ellis-Jaffe sum rules are related to polarized deep inelastic electron-nucleon 
scattering, which is described by the hadronic tensor 

W^, = j-Jd^ze"i^p,s\J^{z)MO)\p,s) 

= [-9,. + '-f) F,ix, Q^) + (p, - {p. - ^g.) Q^) 

+^e^^.paqp ( ^^H^' Q ) + f^T^p 92{x, Q)j ■ (2) 

Here = ^i'lpilp^i is the electromagnetic quark current where Cj = 2/3, —1/3, —1/3, ■ ■ • 
is the electromagnetic charge of a quark with the corresponding flavour, u, d, s. x = Q^/ {2p- 
q) is the Bjorken scaling variable and = —q^ is the square of the transferred momentum. 
\p,s) is the nucleon state that is normalized as {p,s\p',s') — 2p°(27r)^5^^^(p — p')Sss'- The 
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polarization vector of the nucleon is expressed as — U{p,s)^albU{pis) where U(p,s) is 
the nucleon spinor U{p, s)U{p, s) = 2mN. 
The Ellis- Jaffe sum rule then reads 



(3) 



8 



9Q2 




and the Bjorken sum rule is the difference between the proton and neutron sum rule: 



where the definitions are the following: 



4 1 



4 

Ins = yC'^; = 7jv5 + 3^/; 



:2k , , 

2/3o / 1 \ 

(i + oK))+o(^^j , 

(4) 
(5) 



Cx{Q) = 1 + J:C 



X 



a, 



(6) 



s=0 



and X = {B, EJ, GLS}, the latter to be defined below. For the Bjorken sum rule, the 
as correction [14], the correction [15], and the correction [16] have been calculated 
in the leading twist approximation. Higher twist corrections have also been calculated [17]. 
The EUis- Jaffe sum rule for the proton and neutron was calculated to order ag [18], to order 
[19] and to order [20] in the leading twist approximation. Power corrections were 
calculated in [21]. The LO renormalization group running of the twist-4 operators have been 
computed in Ref. [22]. The non-perturbative matrix elements are defined in the following 
way: 

\gA\sa = 2(p, s|J^''^|p, s) = {Au-Ad)s^, 

ass^ = 2V3{p,s\4'^\p,s) = (Au + Ad-2As)s^, (7) 
ao{lJ''^)sa = {p,s\J^\p,s) = (Au + Ad + As)sa. 

where J^'"(x) = '4>^„^^t""4){x) is the non-singlet axial current, where t"- is a generator of the 
flavour group, and J^(.x) = Y^=i V''i7o-75^i(^) is the singlet axial current. \gA\ is the absolute 
value of the constant of the neutron beta-decay, \gA/ 9v\ = F + D = 1.2695 ± 0.0029 [23]. 
as = 3F — D = 0.572 ± 0.019 is the hyperon decay constant^. The matrix element of the 
singlet axial current ao(/i^) will be redefined in a proper invariant way as a constant oq: 



oo = exp 



(8) 



^We will obtain this number from Hyperon decays, see [24]. F/D = 2/3 in the large Nc. 
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We use the notation /S.q{jj?)sa = {p, s|g7o-75(?|p, s), q = u,d, s, for the polarized quark distri- 
butions, /o, fs and /§ are the twist-4 counter parts of Oq, CI3 and Og. /j's are scale dependent 
and here they are defined at Qq, i.e. fi is the reduced matrix element of i?2(T5 renormalized 
at Qg, which is defined for the general flavor indices, with f being the flavor matrices, as 



Ri„ = gipG^^Yt'ip, {p, s\RiJp, s) = fiS^ 



0, 



(9) 



and G^i, = \e^yai3G"^ is the dual field strength. 

In the left-hand side of equations (0111) target-mass effects have been included using the 
Nachtmann variable [25]. They read 

X 1 m^x^ C, 

^ ~ 9 g2 - 



gi{x, Q^)dx + 



o 



/if 



(10) 



where ^ = 2x/ (l + y'l + 4m^a;^/ is the Nachtmann scaling variable, m-Ar is the nucleon 
mass. The quantity Mi is the first Nachtmann moment of gi that absorbs all the target 



mass corrections, 



ijn\/Q 



2\n 



and 



m 



N 
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(af + Ad 



where is the target mass correction given by the x -weighted moment of the leading-twist 
gi structure function, and c?^ is a twist-3 matrix element given by 



d^ 



dx x^ (2c/f + 3^, 



(12) 



We will consider the Ellis- Jaffe proton sum rule and the Bjorken sum rule. For the former 
we will use the data points given in Ref. [26], which are already given in terms of Mf , and 
for the latter we will use the data points given in Ref. [27], which used the values^: 



-0.0029. 



p—n 
a'2 



0.0279 . 



(13) 



If one also considers DIS of neutrinos with nucleons, the GLS sum rule appears, for which 
the leading twist have been computed with NNLO accuracy (see for instance [28]): 



M^^s ^ 1 (Mf (Q2) + M-{Q')) = 3Cgl5(«s 



8 fasiQl 



where 



M3^(Q^ 



9 \as{Q^) 



- C dxFf{x-Q^)^ 







21 + 



4m%x'^ 



dx F^^(x;Q^ 

Jo 



2 mjfX 



^ 3" g2 



+ 



m 



N 



1/2- 



(14) 



(15) 



^To take them as constants is an approximation, nevertheless, their effect on the fit is small in comparison 
with other source of errors. 
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{P\Ol\P).^in ^ 2P,{{Ol)) ^ 2P,fi , (16) 

averaged 

and 

Ol^ = uG^.-f'^j^u + dG^,Yl5d. (17) 

For this sum rule, wc will use the data of the CCFR collaboration [29]. 

The difference between Cgls and Cb first start at 0{al) and is proportional to the 
number of light fermions. This new term is of a "light- by-light" nature and proportional to 
a new casimir. The anomalous dimension of the higher twist is equal to the Bjorken case. 

The difference between the three Cx is the light flavour dependence. In the limit n/ — > 
they are all equal. 

We also consider the correction due to the charm quark (with finite mass) to the pertur- 
bative series of Cb and Cgls- They have been computed in Ref. [30]. The 0{a'^) correction 
is equal for both of them and rather small (actually negligible compared with other source 
of errors). Note however that the leading order correction is different in each case (zero 
for the Bjorken case). This correction depends on the Cabibbo angle. We take the value 
sin 6* = 0.224 [23]. According to Ref. [30], it is a good approximation to work with 3 light 
flavours plus one massive flavour up to rather large energies. This is the situation we will 
consider in this paper. 



3 Renormalons 

The Wilson coefficients Cx can be expressed in terms of Sx, their Borel transform, defined 
as 

°° r^") //Lit \" 

in the following way: 

Cx{Q) = l + -r Sx{u)e^^^''du . (19) 

However, the perturbative expansions in ctg of the Wilson coefficients Cx are expected to be 
asymptotic and non-borel sumable. In other words, we expect to have singularities in the 
real axis of Sx{u). Those lying in the positive axis are called infrared renormalons, and those 
lying on the negative axis are called ultraviolet renormalons. The position and strength of 
the singularities can be obtained by using the renormalization group and consistency with 
the operator product expansion. In particular, the infrared renormalons of the perturba- 
tive series are obtained by demanding their cancellation with the ultraviolet renormalons 
of the higher-twist terms. Although the renormalon cancellation has only been explicitly 
shown in some cases in the large-/3o limit, it is assumed to hold beyond this approximation. 
Based on this assumption, one may obtain additional information on the structure of the 
infrared renormalon singularities of the matching coefficients, based on the knowledge of 
the ultraviolet renormalons in higher- dimensional matrix elements, which are controlled by 
the renormalization group [5]. This model- independent approach has been applied in heavy 
quark effective theory in [10, 11,31,32]. In our case, the singularities closest to the origin are 
located in the real axis at u — ±1. 
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The ultraviolet renormalon structure of the moments of the DIS structure functions have 
been computed in Ref. [33]. For those one gets the ultraviolet renormalon for the sum rules 
we are discussing here, which is the same in all cases up to a constant. For the case of 
the GLS or Bjorken sum rule, the ultraviolet renormalon formally dominates for Uf > 2 for 
n — »• oo. For the Ellis- Jaffe sum rule, since the infrared renormalon is weaker, the ultraviolet 
renormalon dominance appears for even a smaller number of flavours. Nevertheless, at low 
orders in perturbation theory the infrared renormalon appears to be dominant. This can be 
seen from the fact that the sign of the known terms of the perturbative series is equal whereas 
if the ultraviolet renormalon were to be dominant we would find a sign alternating series. 
Nevertheless, we will perform a conformal mapping to avoid the ultraviolet renormalon. The 
fact that we will obtain similar number than without conformal mapping will support the 
view that the normalization constant of the ultraviolet renormalon is small in comparison 
with the infrared one. Indeed, a similar conclusion was obtained in Ref. [34] using Fade 
approximants. Therefore, we will neglect ultraviolet renormalon effects in the leading-twist 
Wilson coefficients in what follows. 

The Borel transform near the closest infrared renormalon singularity has the following 
structure: 

Sxiu) = ^ivj/^)^— — + df{a -u) + d^{a - u)' + ■ ■ ■) + Sre,{u) , (20) 

where Sreg{u) is an analytic function at m = a and we define b = Pi/Pq. The procedure 
to fix the coefficients of this expansion (except Nx^^) is by demanding consistency with 
the operator product expansion. In other words, we demand the ambiguity of the Borel 
transform to cancel with the ambiguity of the ultraviolet renormalons of the twist-4 matrix 
elements (see Eqs. (j3l4ll4|) ). Therefore, 



Im 



oo n J^? 







Sx{u)efn°">''''^^du 



oc ^«,(Q)-^- 



^^el^a^iu)-'-'- (l - ^a^{u) In (^] ) . (21) 



This fixes a = 1 and bx [6]: 



^GLS - Ob - -ITTT ' OEJ - -TTiT 



(22) 



bx dictates the strength of the singularity. It is interesting to study its dependence on Uf. 
In the Bjorken and GLS sum rules, for Uf G (0,6) =^ 1 + b + bx € (1,2) so, formally, one 
could just keep the first two terms of the series in Eq. since the next term would go to 
zero for m — >■ 1. This is also the case for Ellis- Jaffe ii Uf < 4, otherwise one could even stick 
to the first term only. 

If the Wilson coefficients multiplying the higher twist operators were known exactly, we 
could also fix the coefficients . Unfortunately, we only know their leading log running. 
Nevertheless, by performing the matching at a generic scale z/, we will be able to resum the 
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terms of the type (1 — m)" ln"((5^/z^^) and obtain the logarithmically dominant contribution 
to ~ ln"(QVi/2). We obtain 

Sx{u) = ^Arj/^)^^_l^iFi {-bx, -b - bx, (1 - u) HQ^u')) + Sre,{u = 1) . (23) 

The leading asymptotic behaviour of the perturbative series due to the first infrared renor- 
malon reads 



X 



f poY r{n + l + b + bx) ^,x T{l + b + bx) r{l + b + bx + n - s) 



Q^\A7tJ T{l + b + bx) fr'o ' T{l + b + bx-s) T{l + b + bx+n) ' 

(24) 

where (1^=1. The logarithmically enhanced contribution to df is known. By introducing 
it to Eq. (j24|) . we obtain 

^ //JoV Tjl + bx + b + n 
Q2 UnJ T{l + b + bx 



Cx = Nxjy^{-\ w.^,^,^ ^F,[-bx,-b-bx-nM{Q /u)). (25) 



The above expression contains subleading terms in the 1/n expansion. In the strict 1/n 
expansion, it simplifies to: 



This result is numerically not very different from the previous expression for large n. On the 
other hand, it will give more stable results than Eq. ()25|) when working in the RS scheme 
for small n. 

With the above results we can identify the contribution to Cx that comes from the 
renormalon. It reads 



n=n* 



where n* indicates the freedom to add and subtract finite-order contributions in perturbation 
theory. 

At this stage, it is interesting to notice that 5Cx{v) can be written in the following form 

up to subleading terms. This expression will be more convenient for our purposes, since 
it has the same scale dependence on Q as the higher twist contribution. Therefore, it can 
be moved from the leading to the subleading twist term without jeopardizing the Q scale 
dependence predicted by the factorization of scales. 
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3.1 Determination of the normalization constant 



In this subsection we will obtain Nx- Let us momentarily neglect the ultraviolet renormalon. 
If this were the case, we could concentrate on the singularity closest to the origin in the Borel 
plane located at u = 1. Then, we can proceed in analogy with Refs. [11,13,35] and define 
the new function 

oo 

Dx{u) = (1 - uy+'+'-Sxiu) = Y: D^u^ . (29) 

n=0 

One would then obtain Nx from the following identity 

oo 2 

Dx{u = 1) = Y.dP = Nx-^,, (30) 

n=0 ^ 

where the first three terms of the expansion are known. Note that, formally, the expansion 
parameter is m = 1, but one does not really know if there are small factors multiplying the 
powers of u. In practice, the series is quite convergent, and stable numbers for Nx can 
be obtained from most of the sum rules and values of Uf. Nevertheless, we still have the 
problem of the ultraviolet renormalon located at m = — 1. Formally, this renormalon would 
make the series in Eq. ()30|) non-convergent. In order to avoid this problem we will perform 
the conformal mapping [36], 

Hu)= ^ (31) 

This transformation maps the first infrared renormalon to w = 1/3 and all other singularities 
to the unit circle |w| = 1. In the conformal mapping the expansion parameter is w = 1/3. 
In practice the effect of doing the conformal mapping is small, which points to the fact that 
the effect of ultraviolet renormalons is small in comparison with the effect of the infrared 
renormalon located at u = 1. We will only give numbers for the conformal mapping case for 
theoretical reasons. Nevertheless, as we have already mentioned, they will be quite similar 
to the computation without conformal mapping. Our best values for Nx can be found in 
Table ^ They have been computed with NNLO accuracy, after conformal mapping, at the 
scale of minimal sensitivity to the scale variation. The scale dependence of the results as well 
as the convergence is shown in Fig. ^for some selected values of nj. We can see that in most 
cases the scale dependence becomes smother as we go to higher orders. The convergence 
depends on the number of flavours. It is optimal for Uf = 3, which actually happens to 
be the most interesting case from the physical point of view and it deteriorates for large 
Hf. The error quoted in Tabled and represented in Fig. ^by the grey band, stands by the 
maximum between the difference between the NNLO and NLO evaluation at the scale of 
minimal sensitivity of the NNLO evaluation and the difference between the NNLO and NLO 
at the scale of minimal sensitivity for each of them. 

The values of Nb and Ngls are consistent with each other within errors. This is consistent 
with the interpretation that the light-by-light term does not contribute to the renormalon 
as it was done in Ref. [36]. On the other hand our value for A^^ appears to be smaller than 
the number obtained in Ref. [34]. 
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Nb 


Nej 


Ngls 





-0.523± 0.154 


-0.523±0.154 


-0.523±0.154 


1 


-0.487± 0.126 


-0.423±0.120 


-0.479±0.121 


2 


-0.451± 0.101 


-0.291±0.085 


-0.436±0.094 


3 


-0.414± 0.079 


-0.103±0.035 


-0.393±0.070 


4 


-0.378± 0.058 


* 


-0.351±0.059 


5 


-0.343± 0.102 


* 


-0.311±0.134 


6 


-0.311± 0.194 


* 


-0.272±0.232 



Table 1: Values of the infrared renormalon residue of the Bjorken, Ellis- Jaffe and GLS 
leading-twist Wilson coefficient Cx after conformal mapping. " *" means that no stable result 
is obtained. 

We are then able to give some estimates for the coefficients of the perturbative series. 
We provide them in tables |2l 01 IH We should stress that our numbers incorporate the right 
asymptotic behaviour, which is not the case for large-/3o estimates. For the Bjorken and GLS 
sum rules they have been calculated in [37]. The comparison with the exact result works 
reasonable well for rif smaller than 6 for the Bjorken or GLS perturbative series. For large 
Uf the comparison with the exact result gets worse. Note that for rif = 6, the normalization 
constant of the renormalon could be almost compatible with zero if the errors are included. 
This fits with the picture that the renormalon is less important when the number of flavours 
grows and one can reach to the point where the infrared renormalon disappears. For the 
Ellis- Jaffe perturbative series the same discussion applies with nj smaller than 4 (for nf = 4 
no stable determination of the normalization constant can be obtained). Again this ffi with 
the picture that the infrared renormalon becomes weaker when the number of active flavours 
increase. 

We can also compare with other estimates one may find in the literature for the Bjorken 
(GLS) sum rules [34,36,38,39]. We find that our predictions somewhat lie on the upper limit 
of the range of values obtained in these references for C'^b/gls- 

As a final comment it should be mentioned that the introduction of the leading logarithms 
in the renormalon, which is novel, does not actually improve the agreement with the running 
on V of the perturbative coefficients of the perturbative series, when they are known. If this 
is an artifact of the leading-order analysis and would be solved at higher orders remains to 
be seen. 



4 Renormalon subtracted scheme 

In the previous section we have obtained the contribution to the perturbative series of the 
leading-twist Wilson coefficient that produces its leading asymptotic behaviour. This be- 
haviour limits the accuracy that can be obtained for the observable from the perturbative 
series alone, which it will always have an error of the order of the higher twist corrections'^. 

•^Actually, this is so for the order in as for which the difference between the exact and the finite order 
result is minimal. If one goes to higher order in perturbation theory the series will deteriorate and the error 
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(a) Nb, Nej and Ngls with n/ = 0. 
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(c) Nb with n/ ~ 6. 
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(e) Ngls with n/ = 3. 
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/ 
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.4 1.6 

= 3. 
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NLO 
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(f) Nqls with Uf — 6. 



Figure 1: Scale dependence of Nx with conformal mapping for Uf = 0,3,6. The dashed (red) 
line is the LO result, the dotted (blue) line is the NLO result and the continuous (black) line 
the NNLO result. The band represents the error. 
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^/ 





1 


2 


3 


4 


5 


6 




0.523 


0.487 


0.451 


0.414 


0.378 


0.343 


0.311 




0.516 


0.452 


0.391 


0.334 


0.280 


0.228 


0.175 


-cf{Q) 


1.295 


1.045 


0.824 


0.630 


0.461 


0.318 


0.198 


-Cf{Q) 


4.199 


3.149 


2.285 


1.588 


1.042 


0.630 


0.334 


-Cf{Q) 


17.07 


11.94 


8.008 


5.097 


3.023 


1.625 


0.750 


-Cf{Q) 


83.91 


54.77 


34.08 


19.93 


10.75 


5.175 


2.100 



Table 2: Renormalon-hased estimates of the perturbative coefficients for u = Q and for 
different number of flavours. We use the expression from Eq. \2(^) except for for which 
we use the expression from Eq. otherwise the result is 0. 





1 


2 


3 




0.423 


0.291 


0.103 




0.392 


0.250 


0.080 




0.868 


0.478 


0.129 




2.547 


1.253 


0.298 


-c'^kQ) 


9.476 


4.222 


0.896 


-cfAQ) 


42.87 


17.42 


3.336 



Table 3: Renormalon-based estimates of the perturbative coefficients C^j for v = Q and for 
different number of flavours. We do not display the column with Uf = since the numbers 
are equal to the Bjorken case. We use the expression from Eq. \2f)]) except for C^j for which 
we use the expression from Eq. / f^) . otherwise the result is 0. 

Therefore, one can not obtain these higher-twist corrections unless the perturbative series is 
defined with power-like accuracy. The specific value for the higher twist will depend upon 
the specific prescription used. Here we will adapt the procedure used in Refs. [11-13]. In 
those references the non-analytic behaviour in the Borel plane that produced the asymptotic 
behaviour of the perturbative series was subtracted from the perturbative series and added 
to the subleading non-perturbative contributions. Actually the definition of the non-analytic 
piece is ambiguous and analytic terms can always be added. The specific quantity we will 
substract to the perturbative series will be Eq. (j28j) with n* = 1^. This quantity fulfils the 
requirement that its imaginary part cancels the imaginary part of the perturbative series 
and that its dependence on Q complies with the structure of the higher twist contribution. 

^This is equivalent to what was called the RS' scheme in Ref. [11], whereas the case n* — was named 
the RS scheme. In our case here both schemes coincide, since the n — contribution from Eq. (|28|l vanishes 
because we expand the F's in 1/n. Obviously, we could also set n* different from 1 (as far as it is not too 
large). This would be equivalent to a change of scheme. The values of the non-perturbative matrix elements 
would change accordingly. 
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1 


2 


3 


4 


5 


6 




0.479 


0.436 


0.393 


0.351 


0.311 


0.272 




0.449 


0.379 


0.317 


0.260 


0.206 


0.154 




1.029 


0.797 


0.598 


0.428 


0.288 


0.174 




3.100 


2.211 


1.507 


0.967 


0.569 


0.293 




11.75 


7.750 


4.838 


2.807 


1.469 


0.658 




53.92 


32.97 


18.92 


9.980 


4.680 


1.841 



Table 4: Renormalon-based estimates of the perturbative coefficients Cq^^ for u = Q and for 
different number of flavours. We do not display the column with Uf = since the numbers 
are equal to the Bjorken case. We use the expression from Eq. \2(^) except for CqI^^ for 
which we use the expression from Eq. \25^) . otherwise the result is 0. 



For illustration, the Bjorken sum rule would read 



9A 



4 1 



where 



and 
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9 

/3,Rs(i^/) = fsii^f) - ^Q^gA 



^0 
'NS 



2/3o 



;i + o(a3))+o(^ 



s=0 



«s(g) 



'NS 



2/3o 



SCb,rs{Q; Vf) . 



(32) 
(33) 

(34) 



.«s('^/). 

Similar changes would apply to the Ellis- Jaffe and GLS sum rules. In Eq. fl33|) and Eq. 
we have subtracted and added the contributions coming from the first infrared renormalon 
at the scale Vf respectively. Thus, in this scheme, the Wilson coefficient Cb,rs(Q; ^f) is free 
of the first infrared renormalon and the associated n\ behaviour. Therefore, the series is 
expected to converge better. On the other hand, the higher twist, fs^Rsi^f), is free of the 
first ultraviolet renormalon. 

In principle, the above series could be improved by incorporating the running in Uf to any 
order in ttg, which is renormalon free and therefore it could be obtained with good accuracy. 
In this situation it is legitimate to use the principal value (PV) prescription (or any other), 
since the renormalon ambiguity will cancel in the ratio. In this situation 6Cx,rs reads 



-bx 



Air 



b+bx 



where 



Db{-x) = -xe-^ {x^cos(7r6)r(-6) - {-xf [T{-b) - T{-b, -x)]} , 



(35) 
(36) 
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and, 

roc 

T{b,x) = / dtt^-^e-\ (37) 



denotes the incomplete F function. One would then work with the following quantity 

c^B,Rs(Q; ^/) = Cb,rs(Q; ^) + [^c^xlkQ; ^) - ^c^Zl^Q- Uf)] , (38) 

where one would set z/ ~ Q and work order by order in perturbation theory in the first 
term in the right-hand side (where there is no large logs), whereas for the other terms the 
PV prescription is used. This allows to partially resum the dependence on z//. Let us also 
note that the first term in Eq. ()36|) corresponds to (up to the anomalous dimension). 

Therefore it cancels in the difference, SCx'XsiQ'i ^) ~ ^^x'RsiQ'i ^f)i ^^"^ ^^'^ ^e neglected. 
Working with the principal value prescription has produced very good results in heavy quark 
physics, where it was possible to check the running of vj with "experiment" (lattice), see [13]. 
However, in that situation, the running in vj was known with a very high accuracy, since 
there was no anomalous dimension and the running on Vf could be deduced simply from the 
running of a^- In our case we only know the leading order and we expect a less accurate 
result. 

Finally, we would like to mention that our knowledge of the running in v of the higher 
twist terms is much more limited than in heavy quark physics. Here, we only know the 
leading log running. This has consequences in Eq. (jHHj) since, once we expand in as('^); there 
are some subleading logs which are unknown. 



5 Comparison with experiment 

In this section we will compare our theoretical predictions with the experimental data. Our 
aim is to perform a combined global fit for the three sum rules. We will use the experimental 
data for the Bjorken sum rule from [27,40-43] analyzed according to Ref. [27]. Note that 
the elastic contribution has to be included in the experimental numbers in order the sum 
rule to be fully inclusive, i.e. 

F(Q2) = \f,{Q') {F,iQ') + F,{Q')) + Fi„ei.(Q') • (39) 

The empirical parameterization of the elastic form factors was taken from [44]. For large 
momentum this contribution is completely negligible. 

We take the experimental data points for Mf((5^) from Ref. [26], which have used the 
experimental results for the structure functions from [45-52]. 

For the GLS sum rule we will use the experimental data from the CCFR collaboration [29]. 

We first illustrate the problem of convergence of the perturbative series by drawing the 
perturbative series in the MS scheme at different orders in in Fig. El The perturbative 
series has a relative good convergence. However, this convergence deteriorates when we 
approach to low energies. We also see how the perturbative theoretical result diverges from 
the experimental numbers at low energies. As we have already stated, the solution to this 
problem comes from using the RS scheme. We plot again the perturbative series in the RS 
scheme in Fig. IHlfor two values of vj: Vf = 0.8 and 1 GeV. We can see that in both cases 
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1 1.5 2 2.5 3 3.5 

Q 

(c) GLS Sum Rule. 



Figure 2: Leading-twist contribution to the 
sum rules at different orders in perturbation 
theory in the MS scheme with u = Q com- 
pared with the experimental data. The long- 
dashed (green) line is the LO result. The 
dashed (magenta) line is the NLO result. 
The dot-dashed (blue) line is the NNLO re- 
sult. The dotted (red) line is the NNNLO re- 
sult and the continuous black line is the ( es- 
timate) N^LO* . At this order we have used a 
renormalon-based estimate for Cx (Q) from 
Tables\B\^\^ ao = 0.141. 



the convergence of the perturbative series improves. This is specially so around the 1 GeV 
region. For vj = 1 GeV we can also see a qualitative change in the figure with a much better 
agreement with experiment around the 1 GeV region. This alone does not mean much, since 
it only reflects the Vf scale dependence of the pure perturbative piece. Nevertheless, this scale 
dependence is known and can be predicted by perturbation theory. This scale dependence 
cancels with the scale dependence of the higher-twist terms. Therefore, the complete result, 
including the higher twist terms, should be independent of Vf. Thus, consistency demands 
that if we perform the fit to experiment including the higher twist fii^f) for different values 
of z//, in particular for 0.8 and 1 GeV, the results obtained for the fiiyf) should be consistent 
with the result obtained by performing the perturbative running in Vf with respect these two 
values. We have checked that this is so within the errors of our evaluation. This reassures 
the reliability of our fit. This also tells us that it is reasonable to use the operator product 
expansion formulas for Q larger than 0.8 GeV. This is the attitude we will take in this paper 
where we will include (unless otherwise indicated) the experimental points for Q > 0.8 GeV. 
Therefore, the perturbative plots in Fig. |21 with Uf = 1 GeV can be reinterpreted as that a 
piece of the higher twist correction has been included in the pure perturbative term. This 
piece corresponds to the running of Uf from 0.8 to 1 GeV of the higher twist term and can 
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B 0.24 



0.2 
. 18 
0.16 
. 14 


. 1 




2 . 8 
2 . 6 
2 . 4 
2 . 2 




1 1.5 2 2.5 3 3.5 

(b) = 1. 



1 1.5 2 2.5 3 3.5 

Q 

(e) I// = 0.8. 




Figure 3: Leading-twist contribution to the sum rules at different orders in perturbation 
theory in the RS scheme with v = Q for two different values of the subtraction point, Uf, 
compared with the experimental data. The long-dashed (green) line is the LO result. The 
dashed (magenta) line is the NLO result. The dot-dashed (blue) line is the NNLO result. 
The dotted (red) line is the NNNLO result and the continuous black line is the (estimate) 
N'^LO*. At this order we have used a renormalon-based estimate for Cx\Q) from Tables\^ 
00 ao = 0.141. 
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be obtained from the renormalization group. Therefore, in this sense, the change of slope 
observed in Fig. |21 can be interpreted as having its origin in perturbation theory. Note that 
the change of slope, from the experimental point of view, comes from the elastic term. 

In Fig. El higher-twist effects have not been included. Our next aim is to perform a fit 
of do and the subleading twist matrix elements from the available experimental data. We 
refrain from trying to fit as, since the experimental errors appear to be too large. Actually, 
they will be one of the major source of uncertainty of our analysis. We perform a global 
fit of all the available data from the different sum rules at the same time. The size of the 
experimental errors is the largest for the GLS sum rule, whereas the most accurate data 
come from the Mf (Q^) experimental points. In any case, will be obtained from the GLS 
sum rule alone, since this fit is independent of the other two sum rules. We perform the 
global fit to different orders in the expansion in ag of the leading twist perturbative Wilson 
coefficient. We work with the running ag consistent with the accuracy one is working at each 
order^. We show the results of the fit at different orders in perturbation theory in Table El 
We can see how the series shows convergence. We can then obtain relatively good estimates 
for do and /s, for the other non-perturbative parameters the situation is less conclusive. The 
values have been obtained with Vf = 1 GeV. 





do 


/o,Rs (1 GeV) 


/3,Rs (1 GeV) 


(1 GeV) 


/8,Rs (1 GeV) 


LO 


0.043 


-1.571 


0.119 


1.030 


5.90 


NLO 


0.120 


-0.817 


-0.075 


0.304 


3.11 


NNLO 


0.136 


0.014 


-0.106 


0.097 


-0.18 


NNNLO 


0.139 


0.552 


-0.117 


0.014 


-2.34 


N^LO* 


0.141 


0.790 


-0.124 


-0.029 


-3.30 



Table 5: Determination of clq and the higher-twist non-perturbative parameters from the 
global fit to the sum rules at different orders in perturbation theory. 

In order to estimate the errors, we allow for a variation of as{Mz) = 0.118 ± 0.003, of 
Nx (according to the error given in Table Q), of the allowed set of experimental points (we 
consider two situations: a) all the data points for q > 0.8 GeV and b) all the data points for 
g > 1 GeV; our central values will be the ones obtained with option a)), and also consider 
the experimental errors. We also consider the difference between the N^LO* and NNLO 
result as an estimate of the error due to the convergence of the series. This error happens to 
be very small in comparison with the other source of errors. Summarizing, we obtain (with 
the fi in units of GeV^) 

ao = 0-141±ng^(^c. 3)lg:gg^(^iVx)j:ni(exp)^°:°°l(pert)^°:°l°(g > 1) , (40) 

^We have also performed the analysis using the four-loop running as at any order. The final results 
are very similar. The use of the four- loop running as somewhat accelerates/improves the convergence of 
the series. Nevertheless, we prefer to keep ourselves consistent and only resum the logs associated to each 
order in perturbation theory. The fit using the principal value prescription is consistent with the finite order 
computation. Nevertheless, it is less precise because of the reasons mentioned in sec. 01 Therefore, we will 
not consider it further in this analysis. 



17 



M 7Qn+0-24Vj:^ N -0.489 
U-'yU_o.399(,()asj +0.489 



X 



-0.124 



-0.050/ 
}-0.032l 



/o,Rs(l GeV) 
/3,Rs(l GeV) 
/5Ss(l GeV) 
/8,Rs(l GeV) 

If we combine all the errors in quadrature we obtain 



^S:}^^(exp)l°Ji(pert)^i:°^°(g > 1] 



(41) 



-0.029"°-^^^ 



{6a 

+0.124(^'^ 



:s;+Hl^(5iv^);S:i^}(exp)inS?(pert): 



;n^^(5A^x);SJll(exp)l°;°^^(pert) ± 0(g > 1) 



-3.30;n?(5«s)l}|^(5iVx);°|?(exp)l°|^(pert)+lJ?(g > 1) . 



ao 

/o,Rs(l GeV) 
/3,Rs(l GeV) 
fi,Ksa GeV) 
/8,Rs(l GeV) 



0.141 ±0.089, 
0.790l}:i^ GeV^ , 

-0.124l° }^7 Q^y2 

-0.029l°:i^ GeV^ 
-S.SOtlfo GeV^ , 



We have also performed the fit with Uf = 0.8. We obtain in this case 



ao 

/o,Rs(0.8 GeV) 
/3,Rs(0.8 GeV) 
/5^Rs(0.8 GeV) 
/8,Rs(0.8 GeV) 



0.143, 

0.365, 

-0.210^ 

-0.167; 

-1.58, 



(43) 
(44) 



(45) 
(46) 
(47) 
(48) 
(49) 



(50) 
(51) 
(52) 
(53) 
(54) 



whereas the magnitude of the errors is similar to the fit with Uf = 1 GeV. As expected, the 
value of oq almost remain independent of Uf. The higher twist parameters do depend on Uf 
in a way predicted by the renormalization group. For instance, for we would have the 
following expression^ 



«s(^'/) 



'NS 



2/3o 



+ 



9 



(55) 

and analogously for the other higher twist terms. Therefore, we should recover the values 
in Eqs. fl45ll49|l after performing the renormalization group running of Eqs. ()5HI54|1 . If we 
perform such running we obtain 



/o,Rs(l GeV) 
/3,Rs(l GeV) 
/5^Rs(l GeV) 
/8,Rs(l GeV) 



0.335 GeV^ 
-0.130GeV2 
-0.004 GeV^ 
-1.46 GeV^ 



(56) 
(57) 
(58) 
(59) 



We see that the running goes in the right direction for /s and whereas for /o and /8 the 
values remain constant. Either way, the numbers agree with those obtained in Eqs. 
within errors. 



^5Cb,ks has a renormalon ambiguity, which cancels in the difference we find in the second term in the 
right-hand side of the equation. We would like to remind the reader that in order to enforce this renormalon 
cancellation order by order in as both terms have to be expanded with as taken at the very same scale. 
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(a) Bjorken Sum Rule. 




12 3 4 

(b) Ellis- Jaffe Sum Rule. 
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(c) GLS Sum Rule. 
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Figure 4: Global fit to the Sum Rules in the RS scheme with v = Q at different orders in 
perturbation theory including the leading and next-to-leading twist term. The long-dashed 
(green) line is the LO result. The dashed (magenta) line is the NLO result. The dot-dashed 
(blue) line is the NNLO result. The dotted (red) line is the NNNLO result and the continuous 
black line is the (estimate) N'^LO*. At this order we have used a renormalon-based estimate 
for Cx\Q) from Tables\^ [Tj The values of ao, fi and f^ are taken from Eqs. 
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Finally, we have also considered the inclusion of l/Q^ corrections. The values of /o and 
/g are not stable under the inclusion of these effects, although there is a correlation on the 
values of /o and /§: a large positive value of /o is only possible if we also have a large negative 
value of fs- The point is that, even if formally it should be possible to distinguish /o and 
fs due to the different anomalous dimension, in practice we do not have enough accuracy. 
Therefore, the values obtained above for /o and fs (and its errors) should be taken with 
caution. This warning also applies to the determination of , for which the inclusion of 
corrections significantly changes its value. For the other coefficients, Sq and /s the 
variations are smaller than the errors of our fit. 

The error appears to be dominated by the experimental one in do and /s, the objects we 
can compute with better accuracy. 

We would also like to note that, in some cases, the values of the higher- twist non- 
perturbative parameters are compatible with zero within errors. 

To illustrate the quality of the fit we plot our final results, the sum rules including the 
leading and sub leading twist, at different orders in perturbation theory with our best fit, 
Eqs. ()45H49p . compared with the experimental data in Fig. |3] 

The analysis of the GLS sum rule taking into account renormalon effects has been con- 
sidered in Ref. [36]. In this reference only two data points were used and a Principal Value- 
like Borel resummation prescription was used. The authors also included some pure non- 
perturbative effects arguing that they could be inherited from the renormalon computation. 
This is still an open question and, actually, it has been criticised by one of the authors in 
Ref. [53]. 

6 Conclusions 

We have studied the large order behaviour in perturbation theory of the Bjorken, Ellis- 
Jaffe and GLS sum rules. In particular, we have considered their first infrared renormalons, 
for which we have obtained their analytic structure with logarithmic accuracy and also 
an approximate determination of their normalization constant. Estimates of higher order 
terms of the perturbative series are given. The RS scheme has been worked out for these 
observables and compared with experimental data. The convergence of the perturbative 
series greatly improves in this scheme, specially around the 1 GeV region. In particular, 
for i^f = 1 GeV, the agreement between the pure perturbative contribution and experiment 
is quite good. The fact that we have a convergent series in perturbation theory allows to 
give meaningful values for the higher twist condensates with well defined errors. We have 
performed a detailed analysis, being able to give predictions for cio and some higher twist 
condensates, including error bars. Our best fits for the sum rules can be found in Fig. EJ 
The values for the non-perturbative matrix elements read 

ao = 0.141 ±0.089, (60) 
/3,Rs(lGeV) = -0.124+°;}37GeV^ (61) 

The experimental situation is not very good for the GLS sum rule, for which we can not 
give a precise number for the higher twist. The experimental precision is not good enough 
to check the assumption of [54] that and are equal. We also do not display here the 
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values of /o and fs, since they have large errors and their values in the fit are somewhat 
correlated. A large value of one of them could be obtained in the fit to the price of having 
the other being large with opposite sign. 

One of the most important source of error of the present analysis is the experimental 
one. Any improvement in this respect will immediately lead to a reduction of the errors of 
the numbers obtained in this paper. 

The quality of the analysis is worse than the one obtained for heavy quark physics analy- 
sis. The sensitivity to the renormalon is smaller and the determination of the normalization 
constant is less accurate that in that case. This was to be expected since the singularities 
in the Borel plane are more far away here than in heavy quark physics. In any case, the 
inclusion of the renormalon cancellation introduces a qualitative change on the perturbative 
behaviour around the 1 GeV region making it much closer to the experimental figures. 

Another issue we would like to mention is that the resummation of renormalon-related 
logarithms does not appear to improve the convergence of the series. On the other hand, we 
have only performed the leading log resummation in this paper. It would be interesting to 
see what happens at higher orders. 

Finally, the possibility to merge with the chiral limit seems closer now but the gap still 
exists. In particular one should find a systematic way to incorporate higher twist effects. 
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